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ABSTRACT 

Filament eruptions and hard X-ray (HXR) source motions are commonly observed in solar flares, which pro- 
vides critical information on the coronal magnetic reconnection. This Letter reports an event on 2005 January 
15, in which we found an asymmetric filament eruption and a subsequent coronal mass ejection together with 
complicated motions of HXR sources during the GOES-class X2.6 flare. The HXR sources initially converge 
to the magnetic polarity inversion line (PIL), and then move in directions either parallel or perpendicular to 
the PIL depending on the local field configuration. We distinguish the evolution of the HXR source motion in 
four phases and associate each of them with distinct regions of coronal magnetic fields as reconstructed using 
a non-linear force-free field extrapolation. It is found that the magnetic reconnection proceeds along the PIL 
toward the regions where the overlying field decreases with height more rapidly. It is also found that not only 
the perpendicular but the parallel motion of the HXR sources correlates well with the HXR lightcurve. These 
results are discussed in favor of the torus instability as an important factor in the eruptive process. 
Subject headings: Sun: activity — Sun: flares — Sun: coronal mass ejections (CMEs) — Sun: X-rays, gamma 
rays — magnetic fields 
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1. INTRODUCTION 

During solar flares, Ha ribbons form along the magnetic 
polarity inversion line (PIL) and separate f rom each ot her in 
the direction perpendicular to the local PIL (|Zirmlll988[) . This 
well-known behavior has been regarded as a piece of evi- 
dence for the so-called standard solar flare model, in which 
the magnetic reconnection proceeds into the higher corona 
(e.g. JPriest & Forbes 2000). Motions of flare ribbons parallel 
to the PIL are also commonly observed but had not received 
much attention until similar motion s are found in hard X- 
ray (HXR) or EUV/UV emissions (Fletcher & Hudsorfe OOll: 
iGrigis & Ben3l2005tlYang et al.ll2009l) . Since the HXR emis- 
sions are due to high-energy particle precipitating into the 
chromosphere, the parallel motion could also be associated 
with the primary energy release from the corona, but its im- 
plication on solar eruptions and flare energy release is yet to 
be explored (iLee & Garvll2008b . 

Recently, the motions of flare emissions parallel to the PIL 
were related t o the phenomenon called asymmetri c filament 
eruption (e.g., Tri pathi et al.l 120061: iLiu et al.l [2009). During 
asymmetric eruptions, only one end of the filament erupts up- 
ward with the other end anchored, which can lead to a sequen- 
tial magnetic reconnection along the PIL and may provide an 
important clue to understanding the parallel motion of HXR 
sources. 

In this Letter, we investigate the 2005 January 15 X2.6 flare, 
in which we found both the asymmetric eruption and the par- 
allel HXR source motion. We will discuss the implications of 
the HXR footpoint motions on the eruption using a non-linear 
force-free field (NLFFF) extrapolation from the active region. 

2. OBSERVATIONS AND DATA REDUCTION 

The f3S active region NOAA 10720 lies close to the disk 
center (N15°, W05°) when the 2005 January 15 X2.6 flare 
peaked at 23:02 UT in GOES soft X-ray flux. We used high- 
resolution Ha - 0.8 A images with a pixel scale of ~0.6" 
and a temporal cadence of 1-2 minutes obtained with the Big 



Bear Solar Observatory (BBSO) to monitor the evolution of 
the active-region filament and that of the early flare kernels. 
The images of the eruptive filament as well as the surround- 
ing coronal structure were also taken at 195 A (with 5.26" 
resolution and ~12 minutes cadence) by t he EUV Imaging 
Telescope (EIT; Delaboudi niere et aT1ll995|) . 

The evolution of the flare HXR emission was entirely 
registered by Reuven Ram aty Hig h Energy Solar Spectro- 
scopic Imager (RHE SSI; iLin etal. 2002). CLEAN images 
(Hur ford et al.l l2002h in the 40-100 keV energy range were 
reconstructed using the front segments of detectors 3-9 (giv- 
ing an FWHM resolution of ~9.8") with 20 s integration time, 
starting from 22:32:04 UT, except the time interval of ±4 s at 
22:34:08 UT when the RHESSI attenuator switched between 
Al and A3 status. Considering the characteristic shape of 
the X-ray sources, we computed the centroid position of each 
HXR source by fitting a 2D elliptical Gaussian above a mini- 
mum flux value, and we estimated the statistical position error 
by varying the threshold at 40%-60% of the maximum flux. 
The footpoint velocity was then derived using three-point La- 
grangian interpolation implemented in the DERIV procedure 
of IDL, as it is less sensitive to the uncertainties in the position 
measurement. The total flux of each individual source was ob- 
tained by summing up the pixel values of the entire source fea- 
ture, and its uncertainty was evaluated as | of the maximum 
flux outside of the source within the field of view (FOV), as 
implemented in the standard RHESSI software package. 

In this multiwavelength study, we used the full-disk 96 
minute magnetogram at 20:51 UT acquired b y the Michelson 
Doppler Imager (MDI; IScherrer et al.l Il995h for data align- 
ment. The accuracy of feature matching is estimated to be 
about one MDI pixel size of ~2". Vector magnetograms of 
the whole active region were t aken by the Digital V ector Mag- 
netograph system at BBSO (Spiroc ket al.l 12002) and were 
calibrated using MDI level 1.8.2 data. T he vector magne - 
togram data have been processed following I Jing et all ([2009): 
(1) the central umbral fields were filled with those measured 
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Fig. 1. — Time sequence of BBSO Ha blue-wing images showing the 
asymmetric filament eruption. The levels of MDI contours superimposed 
in the top panel are ±800, 1600, and 2200 G. The thick line is the main PIL. 
All the images presented in this paper are aligned with respect to 20:51 UT. 

by MDI to alleviate the influence of polarization saturation; 
(2) the 180° azimuthal ambiguity in the transv erse fields was 
resolv ed using the "minimum energy" method (Metcalf et al. 
2006); (3) the projection effects were removed by transform- 
ing the observed fields to heliographic coordinates; and (4) the 
Lorentz for ces and torques were min imized by a preprocess- 
ing method ( Wiegelmann et al. 2006). The NLFFF model was 
then constructed within a box of 248 x 248 x 248 grid points 
with a s ize element of \" using the weighted op timization 
metho d (IWiegelmannl2004l) adapted to BBSO data dJing et al.1 
2009), the results of which show similar magnetic structure 
during —20:30-21:30 UT. For the purpose of this study, we 
analyzed the preflare NLFFF model at 20:51 UT, which was 
chosen for the optimal observing condition. 

3. EVENT EVOLUTION 

Based on the HXR lightcurve, we divide the whole event 
into five phases, I-V (see Fig. [3)2). We describe characteristic 
flare activities in each phase in the following subsections. 

3.1. Asymmetric Eruption (Phase I) 

Time sequence of Ha blue- wing images in Figure Q] shows 
that the western part of the filament (/) lying along the PIL 
gradually rises upward from west to east during —21:41- 
22:00 UT and subsequently erupts, while its eastern part re- 
mains undisturbed as can be seen in the time-lapse movie. 
This filament eruption from the western edge of the ac- 
tive region produced a M1.0 flare starting at 22:01 UT. The 
M1.0 flare appears to be a confined flare, since its flare rib- 
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Fig. 2. — Time sequence of SOHO/EYT 195 A images showing the erup- 
tions of the filament / and the subsequent CME. The dashed box in the upper 
left panel indicates the FOV of the bottom 8 panels of Fig.[T] 

bons/kernels remain almost fixed till the occurrence of the 
subsequent X2.6 flare at 22:24 UT (cf. images at 22:04 and 
22:24 UT). The c onfined flaring a lso implies that / undergoes 
a failed eruption ( Schrijver 2009). 

In Figure [3 we further examine the following eruption 
based on the evolution of overlying coronal loops shown in 
EUV images. Although / erupts (see images at 22:00 UT), 
most probably it does not immediately open the overlying ar- 
cade, which expands from 22:00 to 22:24 UT (see the differ- 
ence image at 22:24 UT) and finally erupts at —22: 36 UT to 
become a fast halo CME (also see iLiu et~aLll2010|) . Hence 
we believe that although the first eruption of / (—21:41- 
22:00 UT) failed, the later eruption (> 22:00 UT) was suc- 
cessful to completely tear open the overlying field, resulting 
in a fast CME and the X2.6 flare. This pr ocess resembles the 
asymmetric filament eruption reported by lLiu et all (2009). 

3.2. Motion of HXR Sources of the X2. 6 Flare 

We show, in Figure the evolution of the HXR footpoint 
sources throughout the X2.6 flare superposed on the Ha blue- 
wing images at the nearest time and the preflare MDI magne- 
togram. We further characterize the motion of HXR sources 
in Figure [4] after a detailed quantitative analysis of the HXR 
sources. 

3.2.1. Converging Motion ( Phase II) 

HXR images can be reconstructed only after —22:32 UT. 
At the event onset during -22:25-22:33:34 UT, the flare ker- 
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Fig. 3. — RHESSI HXR light curve and images together with their centroid positions (pluses) showing the event evolution during four distinct phases (II-V). 
Contour levels are 30%, 50%, 70%, and 9 0% o f the maximum flux at each instance. The black bordered (b) and dotted (c-f) lines are the main PIL and our 
tracing of the filament / (see Fig. [I}. See § !3.2l for details. 

nels Nl and SI show a consistent converging motion in Ha 
blue wing and HXRs (Figs.0?c and|4]i). Meanwhile, the flare 
shear, defined as the angle between the li ne connecting flare 
kernels and that perpendicular to the PIL dJi et al.ll2007h . de- 
creases steadily (Fig.[4]b). 

Similar behavior was observed before where flare kernels 
converge mainly in the direction perpendicular to the PIL for 
the e ntire rising phase of HXR emission (e.g.,|JLetalj[2006, 
120071) . In contrast, the converging motion in the present flare 
is mostly parallel to the PIL, which is oriented mainly in the 
west-east direction (Fig. [2?). 

3.2.2. Parallel Motion (Phase III) 

During 22:33:34-22:45:02 UT, Nl and SI keep a nearly 
constant separation and both move preferentially parallel to 
the PIL (Figs.[3]W and[4]:<i). The angle of flare shear is seen 
to fluctuate until ~22:38 UT but decreases overall at a higher 
rate than that of phase II (Fig.|4t?). Interestingly, Nl first starts 
at a slower speed (~9.7 km s -1 ) and later the speed increases 
substantially (~56 km s -1 ) as approaching the main HXR 
peak PI (Fig. [4]:). The maximum speed reaches >100 km s -1 
(Fig. [4]?), which is comparable to other reported values of flare 
kernel motions found using high-cade nce Ha blue- wing filter- 
grams (Oiu et al. 2002|;|LeeeLal. 2006). 

In Figure |4^ and 0]/*, we plot the speed u and the quan- 
tity uB 2 of Nl in comparison with the HXR flux of Nl (Jhxr 
in Fig. We focus on Nl, because it has the strongest 
HXR emission and exhibits the most intriguing pattern of mo- 
tion. The quantity uB 2 is regarded as a proxy of magnetic 
energy rel ease rate under th e standar d 2D model (fir st pre- 
sented by llsobe et al.1 120021: see also iLee et al.1 120061) . We 
compute B by averaging the four nearest pixels in the pre- 
flare MDI magnetogram around the position of Nl centroid, 
and we present two vector components of u along the east- 
west and the south-north directions as uew and usn, respec- 
tively. Remarkably, five peaks (dotted lines in Fig. [4]) of Ihxr 
have corresponding peaks in uewB 2 , contributed by the foot- 
point motion parallel to the PIL (uew), while only one peak 
of Ihxr at ^22:38:22 UT (dashed line in Fig. [4]) is cotempo- 
ral with that of usnB 2 due to the motion perpendicular to the 
PIL (usn)- This peak-to-peak correspondence is significant 
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Start Time (15-Jan-05 22:24:00) 

Fig. 4. — Temporal evolution of derived parameters for the HXR sources 
(phase III shaded in grey), (a) HXR flux and photospheric magnetic field 
strength at the source centroid of Nl obtained by MDI with a largest noise 
level of 30 G. (b) Angle of flare shear, (c) East-west and south-north motions 
of Nl. (d) Separation of sources Nl and SI. (e) Velocity of motion of Nl. 
(f) Product uB 2 . Vertical dotted and dashed lines mark the selected peaks of 
HXR flux. 
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Fig. 5. — NLFFF model computed using a preflare vector magnetogram (a), and the derived distribution of the decay index along the PIL/filament (b) and in 
different flare phases (c). See §2]for details. 



as it remains when a three-point smoothing is applied on the 
displacement profile prior to the derivation of the footpoint 
velocity. In phase III, the calculated Pearson correlation coef- 
ficient r indicates a strong correlation between I H xr and uew 
with r = 0.71 and a weak correlation between Ihxr and uewB 2 
with r = 0.44, but only poor correlation is present between 
Ihxr and usn or usnB 2 . Similarly for SI, the HXR flux bet- 
ter correlates with the parallel motion than the perpendicular 
motion. 

As a comparison, Gri gis & Benzl (|2005|) also reported a mo- 
tion of two HXR footpoints parallel to the PIL. In their event, 
the parallel motion was rather smooth and no correlation be- 
tween the motion and the HXR flux was found. The au- 
thors proposed that the motion of HXR sources simply re- 
flects 3_movmg, trigger propagating continuously along the 
PIL. iLee & Garyl {2008) presented a generalized framework 
in which both the parallel and perpendicular motions rep- 
resent the magnetic flux change in the corona as a way to 
explain the correlation between the flux of a coronal HXR 
source and parallel motions of Ha ribbons. In the present 
case, the parallel motion of HXR footpoint sources shows a 
good correlation with the HXR flux. This suggests that the 
parallel motion not only maps the propagation of the trigger 
but also represents the primary flare energy release. 

3.2.3. Mixed Motion (Phase IV and V) 

Phase IV includes the peak PI of HXR emission (22:45:02- 
22:57:42 UT), during which Nl and SI separate from each 
other (Figs. \5j?e and|4]i) and Nl turns northward perpendicu- 
lar to the PIL (Fig. |4]:). Ihxr has a main contribution from the 
south-north motion (usn) and shows a moderate overall corre- 
lation with usnB 2 with r = 0.5, in agreement with the standard 
model. The flare shear first decreases rapidly during the HXR 
maximum before 22:50:40 UT then slows down (Fig.[4]b), cor- 
responding to higher and lower northward speed of Nl, re- 
spectively. However, S 1 does not move systematically at this 
ph ase. Similar type of m otion was found in an event studied 
by iKrucker et al.l (|2003 ) in which one HXR footpoint source 
moves roughly along the PIL and is generally correlated with 
HXR flux, while other two footpoints do not show any sys- 
tematic motion. 



Phase V has two minor HXR peaks P2 and P3 (after 
22:57:42 UT), with the conjugate footpoints, N2 and S2, 
showing up in distinct, eastern section of the PIL (Fig. Et/). 
The motion of S2 is parallel to the PIL, while that of N2 shows 
both parallel and perpendicular components. A general cor- 
relation exists between all these motions with HXR flux. It is 
also noteworthy that the flare shear actually increases ^30° 
from P2 to P3. This may be linked to har der HXR spectra of 
P3 than P2 found by Saldan ha et al.l J2008) for the same event. 

4. MAGNETIC FIELD STRUCTURE 

For the eruption mechanism, we consider the torus instabil- 
ity (TI) the most appropriate for this event, because there was 
a continuous flux emergence in the region (IZhao et al.l l2QQ8) 
and the supply of twiste d magnetic flux into the corona tends 
to cause the TI (e.g.. IFan & Gibson! 12007b . Kink instabil- 
ity (iTorok & K liem 2005)) is excluded here because writhing 
motion of the filament is not apparent. Other well-known 
models for solar eruptions , for instance, magnetic break- 
out mo del dAn tiochos 1998) and tether-cutting reconnection 
model (iMoore et al.ll2001|) are not considered because there 
are no obvious chromospheric brightenings before the erup- 
tion, either in remote regions or within the initial flaring site 
at the western edge of the active region. 

In the TI model, the likelihood of eruption is gauged by 
the decay index, defined by n = -d\og(B) / dlogih). Here B is 
the strength of the external field confining the erupting core 
field, h is the height above the photos phere, and n > 1.5-2. 
is the theoretical instability condition ( Kli em & Tor ok 2006). 
We calculate the decay index at 42-105 Mm using B(h) of 
the NLFFF model (§ [2]) as an approximation of the external 
field. Filament activation usually occurs i n the simi lar height 
range according to limb observations (see Liu 2008, and ref- 
erences therein). We note that the NLFFF model used here is 
constructed at a single time in the preflare stage. Ideally we 
should have used time dependent NLFFF models to exam- 
ine the evolving magnetic fields. To obtain reliable extrapo- 
lation, however, we cannot use flare-time magnetograms. We 
therefore only point out which field lines in the pre-emption 
magnetic configuration are involved with the eruption in each 
phase. Figure 02 shows the magnetic field lines stemming 
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from Nl and N2 at their locations in each phase, marked by 
different colors. We calculate n at each position along the 
PIL with index X PIL (see Fig.[5]2 and also Fig. [T] to/?) and plot 
this spatial distribution of n in Figure [3?. At each time of 
Ha/HXR images, we also read off, from Figure 0b, the decay 
index at the particular X PIL where the line joining two conju- 
gate flare footpoints intersects the PIL (see Fig. [3?). In this 
way we are able to approximate the decay index for the coro- 
nal region undergoing magnetic reconnection as a function of 
time, as shown in Figure [5]:. 

In phase I, Ha and EUV images show that the filament / at 
the western edge of the active region failed to erupt, and that 
the associated Ml.O flare is a confined flare. At the location 
above /, the decay index n has an average value of ~1.84. 
Since this value is relatively low, the failed eruption may be 
expected. Although it failed, the Ml.O flare must have af- 
fected the local magnetic field in a way to reduce the overlying 
magnetic flux to some extent, so that the next CME eruption 
can occur along with a stronger X2.6 flare. The successful 
CME eruption is anticipated, as the main flare phases II-V 
progresses from regions of n about 2.5 to 3 with an average 
value of ~2.7. 

In phase II, low-lying and highly- sheared field lines (red 
in Fig. 02) are involved in the reconnection process, during 
which the flare kernels Nl and SI converge with decreasing 
flare shear. This is consistent with the previous view that con- 
verging motion of flare kernels is due to re connection pr o- 
ceeding from more to less sheared field lines (|Ji et al .1 12007). 

In phase III, Nl and S 1 move parallel to the PIL from east to 
west, which corresponds to the motion from strongly sheared 
to less sheared fields (green) and also from stronger to weaker 
confinement of overlying field (Fig. 0:). The location and 
direction of this HXR source motion thus support the idea 
that the CME, which abruptly erupts outward at the same time 
dLiu et al.ll2010|) . acts as the moving trigger to successively 
open the overlying field, in a way similar to the asymmetric 
filament eruption. 

In Phase IV (blue), one footpoint N2 shows a standard per- 
pendicular motion while its conjugate footpoint S2 does not. 
In phase V, N2 and S2 generally move in the opposite direc- 
tion from west to east (yellow). Overall, phase V occurs in the 
region with much enhanced decay index compared to phases 
II-IV. In all phases, the reconnection evolves in the direction 
of increasing decay index, namely, toward the region of the 
maximum n « 3. It thus appears that the magnetic reconnec- 
tion tends to proceed toward the regions with weaker mag- 
netic confinement. 

5. SUMMARY 

We have presented a near disk-center event that shows 
asymmetrical filament/CME eruptions and associated flare 
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kernel motions observed at Ha and HXRs, and discussed the 
results with the aid of NLFFF model. Major results are as 
follows: 

1 . The first activities occurred in the region of low decay in- 
dex implying a strong magnetic confinement, which has pro- 
duced a failed filament eruption and a relatively small flare of 
Ml.O. 

2. The subsequent main event was initiated at the same re- 
gion but proceeded to nearby regions of higher decay index, 
i.e., weaker magnetic confinement. The eruption at this stage 
produced a CME and a stronger flare of X2.6 with compli- 
cated HXR footpoint motions, indicative of complex coronal 
magnetic structure. 

3. Converging HXR footpoint motion was found in the ini- 
tial stage of the X2.6 flare. As suggested before, it is consid- 
ered as a result of reconnection proceeding from more to less 
sheared fields. In contrast to other events, the converging mo- 
tion at the ons et of this flare has a major component parallel 
to the PIL (cf. Ui et al.ll2006l 12007b . and the flaring scenari o 
does not involve quadrupolar reconnection fcf. lJi et aI1 l2008). 

4. During the rising phase of HXRs, footpoints moved 
parallel to the PIL, a phenome n on similar to what w as re- 
ported by iKrucker et alJ (l2003|), iGri gis & Benz (2005), and 
iLee & Garvl (120081) . and bv lYang et alJ (120091) in a statistical 
analysis including this event. Grigis & Benz (2005) suggested 
that the parallel HXR motion occurs because the trigger of 
magnetic reconnection moves along the PIL. In this case, we 
identified the moving trigger with the asymmetric eruption of 
the CME. It is further notable that the parallel HXR source 
motion has a peak-to-peak correlation with the HXR flux in 
complement to the standard 2D model. 

In summary, the observations and interpretations presented 
in this Letter support the idea that an asymmetric eruption 
progressively opens the overlying field from one end of the 
PIL to the other end, which results in a magnetic reconnection 
proceeding along the PIL. A new intriguing property is that 
the asymmetric eruption and thus the corresponding trigger 
of magnetic reconnection tend to progress toward the region 
of weaker magnetic confinement. 
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